home *** CD-ROM | disk | FTP | other *** search
/ MacFormat 1994 August / August CD.bin / Shareware / Education / numericalmethods Folder / chap_10 / trisys.m < prev   
Encoding:
Text File  |  1994-06-05  |  523 b   |  21 lines  |  [MATF/MATL]
  1. function X = trisys(A,D,C,B)
  2. % X = trisys(A,D,C,B)
  3. % Solution of a triangular linear system AX = B.
  4. % It is assumed that D and B have dimension n,
  5. % and that A and C have dimension n-1;
  6. % A sub diagonal, input.
  7. % D diagonal vector, input.
  8. % C super diagonal, input.
  9. % B right hand side vector, input.
  10. % X  solution vector, output.
  11. n = length(B);
  12. for k = 2:n,
  13.   mult = A(k-1)/D(k-1);
  14.   D(k) = D(k) - mult*C(k-1);
  15.   B(k) = B(k) - mult*B(k-1);
  16. end
  17. X(n) = B(n)/D(n);
  18. for k = (n-1):-1:1,
  19.   X(k) = (B(k) - C(k)*X(k+1))/D(k);
  20. end
  21.